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ABSTRACT 



Context. In recent years, the solar chemical abundances have been studied in considerable detail because of discrepant values of 
solar metallicity inferred from different indicators, i.e., on the one hand, the "sub-solar" photospheric abundances resulting from 
J spectroscopic chemical composition analyses with the aid of 3D hydrodynamical models of the solar atmosphere, and, on the other 

Pm , hand, the high metallicity inferred by helioseismology. 

^0 ' Aims. After investigating the solar oxygen abundance using a CO^BOLD 3D hydrodynamical solar model in previous work, we 

undertake a similar approach studying the solar abundance of nitrogen, since this element accounts for a significant fraction of the 
overall solar metallicity, Z. 

Methods. We used a selection of atomic spectral lines to determine the solar nitrogen abundance, relying mainly on equivalent width 
^ . measurements in the literature. We investigate the influence on the abundance analysis, of both deviations from local thermodynamic 

^—^ • equilibrium ("NLTE effects") and photospheric inhomogeneities ("granulation effects"). 

' Results. We recommend use of a solar nitrogen abundance of A(N)=7.86 + 0.12*, whose error bar reflects the line-to-line scatter. 

, Conclusions. The solar metallicity implied by the CO'BOLD-based nitrogen and oxygen abundances is in the range 0.0145 < Z < 

0.0167. This result is a step towards reconciling photospheric abundances with helioseismic constraints on Z. Our most suitable 
estimates are Z = 0.0156 and Z/X = 0.0213. 

^ ' Key words. Sun: abundances - Stars: abundances - Hydrodynamics - Line: formation 

o 

' 1. Introduction C/O ratio cannot be inferred from meteorites; although, isotopic 

■ ratios can be measured with high accuracy. Similar problems af- 

The knowledge of the precise nitrogen abundance in the solar f^^^ nitrogen in an even more severe manner. As a consequence, 

O photosphere is important because N accounts for about 6% (by j^e meteoritic C, N, and O abundances are lower than the so- 

^ mass) of the solai- metallicity, and hence is one of the most im- j^. photospheric ones by large factors; 2-3 for oxygen, about 

O portant contributors to Z after O (53%), C (23%), and Ne (14%). ^j. carbo n, and more than 40 for nitrogen (see e.g.. Palme 

_ Amongst others, the value of Z has a crucial impact on stel- ^ j^^g^^ 2003^). Analyses of observations of the solar corona 

lar structure calculations. The abundance of C, N, and O in ^j^^ pjoyj^e reliable abundance measurements of these 

the solar system can be determined only from the solar spec- elements. Several processes affect the solar corona composition, 

H ti-um. Meteoritic samples are un-representative, because C, N, partially changing the abundance of some elements (Anders & 

and O are volatile elements. With condensation temperatures of Gi-evesse1T989l ; iGrevesse & Sm wal. 2000). The solar niti-ogen 

only a few hundred Kelvin, they condense only partially so that abundance can be derived reliably only from the photospheric 

their abundances in pristine meteorites (e.g., carbonaceous chon- spectrum 
drites) do not provide valid estimates of the average solar system 

composition. The main molecule formed by C and O, cai'bon Besides its importance in determining the solar metalhcity, 
monoxide (CO), is also extremely stable. As a consequence, all niti-ogen is an interesting element from a nucleosynthetic point 
C will be locked into CO if the original material is O-rich, and of view. Its main production channel is the CNO cycle, in which 
will be unable to form solid compounds. Only the remaining C it is produced at the expense of C and O. The main astrophysical 
may condense into grains, and in the resulting meteoritic sample N production sites remain unclear; AGB stars are good candi- 
carbon will be underrepresented. The reverse is true if the envi- dates, in addition to rotating massive stars. Niti-ogen is measured 
ronment is C-rich; all oxygen remains in the gas phase and the in many different objects; hot stars, cool stars, Galactic and ex- 
remaining C may form grains. Therefore, the global elemental tragalactic Hii regions, and Damped Lyo' systems at high red- 
shift. For all of these studies, a good solar reference value is of 



>< 



A(N)= log ;j(N) - log n(H) + 12 fundamental importance. 



The measurements of heavy element abundances inferred 
from photospheric solar spectra should agree with results from 
helioseismology. Helioseismology can provide a measurement 
of the solar metallicity in essentially three ways: i) from the 
depth of the convection zone: this depends sensitively on the 
opacity at the base of the convection zone, wh ich in turn is a 
function of the abundance of heavy elements (iBasu & AntiaL 
11997 ): ii) using information from the core: the small-frequency 
spacings of low-degree modes and their separation ratios are sen- 
sitive to the mean molecular weight in the core, whic h can be re- 
lated to the metallicity of the outer layers (iBasu et al.„ 2007): iii) 
from the sound speed gradient in the ionisation zone: the depth 
profile of this quantity can be inferred from helioseismic inver- 
sions, and comparison with the results of theoretical solar struc- 
ture models o f different met allicity allows us to place tight con- 
straints on Z dAntia & Basuii2006) . A comprehen sive review on 
heliose ismology and solar abundances is given bv lBasu & Antial 
(12008b which the reader is referred to for further details. It is re- 
markable that all helioseismic methods provide results that are 
fairly consistent with each other, which infer a higher metallic- 
ity than currently deduced from the analysis of the photospheric 
spectrum. 

In this paper, we determine the photospheric nitrogen abun- 
dance by using data for atomic lines and the latest CO^BOLD 
solar model. One of our aims is to combine this result with our 
other solar abundance determinations to obtain a measure of the 
solar Z based on the CO^BOLD solar model, and to compare this 
with the Z values inferred by helioseismology. 




Fig. 1. The temporal and horizontal average of the temperature 
profile of the 3D model (solid line) and the temperature profile 
of IDlhd models of three values of aMLT, are shown as a func- 
tion of (Rosseland) optical depth. In addition, equivalent width 
contribution functions (lower, roughly Gaussian-shaped curves) 
for the 746.8 nm line at disc-centre are plotted as a function of 
monochromatic continuum optical depth. 



2. Nitrogen abundance indicators 

2.1. Forbidden lines 

The strongest forbidden [Ni] lines within the ground configu- 
rati on are close to 10 40. Onm and belong to the multiplet ^D°- 
^po. iHouziauxl (Il96lh claimed to detect these lines in the solar 
spectrum, althou gh studies by iLambert & Swings! (1 19671) and 
ISwensson (' 1967) could not confirm these detections. ISwenssonI 
(1967.) also demonstrated convincingly that the predicted equiv- 
alent widths are too weak to be detectable in the Jungfraujoch 
Atlas. 

2.2. Permitted lines 

The permitted N i lines of low excitation potential occur at UV 
wavelengths, mostly in the 80-120nm range. The atomic Ni 
lines measurable in the solar spectrum are found in the visual and 
near IR range, and all have a high excitation potential (Eiow > 
10.3 eV), implying that these lines are formed in the deep layers 
of the photosphere. They are therefor e hardly affected by dep ar- 
tures from LTE (A A(N) ~ -0.05 dex. iRentzsch-HolmL 1 1 9961 ) . In 
the classical analysis that is based on ID model atmospheres, 
these lines have the disadvantage of being very sensitive to the 
temperature structure of the deep layers, which, for theoretical 
models, depends on the treatment of convection and, most sig- 
nificantly, on the choice of the mixing-length parameter (see 
Fig ID- In principle, these difficulties are overcome by mod- 
ern 3D hydrodynamical simulations, which provide a physically 
consistent description of convection from first principles. 

2.3. Molecular lines 



Grevesse et al. (119901) illustrated that atomic and molecular lines 
react in opposite directions to temperature perturbations, so that 



the agreement between the atomic and molecular lines provides 
support for the choice of solar model. In their molecular anal- 
ysis, they and other authors have relied on the use of NH and 
CN molecular lines. In this paper, we derive the solar N abun- 
dance from only the atomic lines and defer the study of molec- 
ular lines to the future. On the one hand the abundance of car- 
bon has not yet been estimated with a CO^BOLD model, and 
on the other ha nd, synthesising the NH A-X band, used e.g. by 
iLamberj (1 19781) in the solar spectrum, would require accounting 
for a large number of lines, both NH lines and blending atomic 
lines, and the state-of the art version of LinforSD is unable to 
handle hundreds of lines. Potential deviations from LTE for the 
molecular species is another issue that hampers the use of molec- 
ular lines for accurate abundance analysis. For these reasons, we 
postpone the study of molecular lines harbouring N. 



2.4. Line selection 

The first selection of lines was completed by extracting all 
the Ni hues fro r ti the NIST Atomic Spectra Database (NIST 
ASD iRalchenkoL 120051) in the range 400-1 200 nm. These lines 
were then checked by visual inspection in the solar atlases of 
Delbouille et al. (1973) and Delbouille et al. (1981); the identifi- 
cation s were achieved b y using the Utrecht spectrum identifica- 
tions ("Moore et all \l96Q for the 400-877 nm range, and those 
bylSwensson et al.l(ll970h for the 900-1 200 nm range. According 
to the NIST ASD, no N i line is present in the 877-900 nm range, 
which is not covered by these identification lists. We also per- 
formed a check of the relativ e line intensities with the Moore 
Multiplet Table (iMooreL 1 19451) . The sti'ongest and least blended 
lines are given in Table[T|with the log gf values, and their qual- 
ities are taken from the NIST ASD. 



Table 1. Selection of permitted N i lines in the visual and near 
IR spectral range 
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Notes: Error of log gf < 0.03 dex (Q=B+), < 0.08 dex (Q=B); niow 
and «up identify the level of the nitrogen model atom used for the NLTE 
computations (cf. Fig.[3ll. 
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Fig. 2. Observed profile of the N i line at A 821.6 nm as extracted 
from the Delbouille (solid) and Neckel (dashed) solar disc-centre 
intensity spectral atlases. The reason for the considerable differ- 
ences in unknown. 



3. Models and line-formation codes 

3.1. 3D hydrodynamical model atmosphere 

Our analysis is based mainly on a 3-dimensional radiation- 
hydrodynamics simula tion (hereafter 3D model) co mputed with 
the CO^BOLD code (iFrevtag et al.', '20021, l2003l; Wedemeyer 
et al.. l2004l) . Some basic information about the setup of this nu- 
merical simulation can be found in iCaffau et al., (2008). who used 
the same model to determine the solar oxygen abundance. A full 
description of the model, including a critical review of its per- 
formance i n reproducing variou s observational constraints, will 
be given in lLudwig et al.l (l2009l) . 

3.2. 1 D reference models 

For comparison, we also used several ID solar models: 

1. The semi-empirical Holweger-Miiller model (iHolweged 
Il967t: Holweger & Muller 1974, hereafter HM); when neces- 
sary, this was placed on a mass column-density scale, using 
the opacity used to produce the CO^BOLD model. 

2. A ID model computed with the LHD code, which uses the 
same microphysics as CO^BOLD and treats convection with 
the mixing-length appro ximation, adopting aMLX^l-O (see 
ICaffau & Ludwig||2007l for further details), and two further 
IDlhd models with cmlt of 0.5 and 2.0, to investigate the 
dependence of the abundance on (Tmlt- 

3. The ATLAS9 solar model with the abunda nces of Asplund 
et al. (2005j as computed by F. Castellfl. 

4. A ID model obtained by temporal and horizontal averaging 
over surfaces of equal (Rosseland) optical depth and over all 
snapshots. We refer to this averaged model as (3D). 

3.3. Line-formation calculations 

The 3D spectrum synthesis computations are all performed 
with Linfor3C0, which can also compute line formation using 
as input different type of ID models. For comparison in the 

' http://wwwuser.oats.inaf.it/castelli/sun/ap00t5777g44377klasp.dat 
^ http://www.aip.de/~mst/Linfor3D/linfor_3D_manual.pdf 



case o f ID mo dels, we also used the SYNTHE code ([Kurucz , 
1993b', '2005b) in its Linux version dSbordone et al.L 12004 : 
Sbordone, 2005) in calculating synthetic spectra. The advantage 
of SYNTHE over LinforSD is that it can easily treat hundreds 
of thousands of lines, while the present version of LinforSD is 
limited to a few tens of lines. 

4. Observational data 

Investigations of solar abundances found in the literature usu- 
ally provide detailed discussions about both the accuracy of the 
adopted model atmospheres and the related uncertainties in the 
derived abundances. However, the accuracy of the observational 
data is often not taken into consideration, since, typically, solar 
abundance determinations rely on a single observed spectrum. 

4. 1 . Differences between high quality solar spectra 

In our study of the solar oxygen abundance (ICafFau et alil2008|) . 
we noted that different high quality solar spectra exhibit differ- 
ences that were larger than expected from their high S/N ratios. 
To investigate the problem related to the differences between the 
solar atlases in the case of nitrogen, we checked four selected 
lines. From this non-exhaustive analysis, we conclude, as we did 
for other elements, that for some, but not all, lines, the solar at- 
lases disagree. This is the case for one of four nitrogen lines that 
we considered as illustrated in Fig.|2] The presence of telluric 
absorption might explain this difference. However, by analysing 
spectra of rapidly rotating hot stars, there is no evidence of any 
telluric contamination, but the day-time sky can exhibit absorp- 
tion that is undetectable at night-time. An indisputable reason 
for the differences still needs to be identified. We strongly be- 
lieve that the astronomical community needs a new high-quality 
solar atlas. 

4.2. Comparison of different EW measurements 

Inspection of the literature often reveals that even the measure- 
ment of equivalent widths (EWs) from one and the same ob- 
served spectrum is less straightforward than what could be ex- 
pected. Different authors use different approaches and measure 



different equivalent widths. A recent example i s the compari- 
son of our KW s for o xygen (ICaffau et al.Ll2008h with those of 
lAsplund et al.l ( |2004) . A long-running debate about the solar 
iron abundance h as in part been driven by differences in the EWs 
(see for example [Holweger et all [1995 , and references therein). 
For the nitrogen, the EW measurements differ also between dif- 
ferent works. Comparing t he measurements of th e EWs of two 
auth oritative works such as Grevesse et al.l (Il990[) and Biemont 
et al. (Il990|) , we conclude that there is good agreement for our 
selected lines of nitrogen. The authors measure similar values 
without any systematic differences being evident. 

4.3. Data used in this study 



In our analysis, we u sed the EWs of iGrevesse eFaP (119901) and 
iBiemont et ai](ll990h , which are in excellent agreement. 

For a few selected lines, we derived in addition the nitro- 
gen abundance from detailed line profile fitting (see Fig.|5l). For 
this purpose, we used t he two centre-dis c intensity atlas es (the 
"Delb ouille" atlas, i.e., iDelbouille eran iT973. Delbouill e et al.1 
ll98lL and the "Neckel" intensity atlas. iNeckel & Labsii 19841) . as 
well as the t wo solar flux atlases (the "Kurucz" solar flux atlas 
Kurucz 2005a and the "Neckel" solar flux atlas iNeckel & Labsl 
1984) ■ Further details about t hese atlases can be found in the 



aforementioned references an d lCaff^au et al] (l2008h . 



5. NLTE computations 

The nitrogen lines considered in our analysis are weak and form 
deep inside the solar photosphere (see Fig.[Tl showing the EW 
contribution function for one of the nitrogen lines). From these 
facts, we do not expect that these lines are very sensitive to 
departures from Local Thermodynamic Equilibrium (LTE). In 
her analysis of t he statistical equilibrium of nitrogen in the Sun, 
iRentzsch-HoimI (17996) inferred that the deviations from LTE are 
small and negative, on average about -0.05 dex. In this paper, we 
confirm this result. 



5.1. 1 D NLTE corrections 

Since the nitrogen lines form at the same photospheric depth 
range in both 3D and ID ((3D), IDlhd) models (see contribu- 
tion functions in Fig.[Tll, this infers that 3D granulation effects 
do not play a fundamental role. Since the contribution functions 
for 3D and ID models are similar, it follows that the computed 
EWs are also similar, implying that the 3D corrections are small. 
Without a code capable of solving the full 3D-NLTE problem for 
nitrogen, we resort to ID-NLTE calculations. 

We computed the departure s from LTE for the (3D) and HM 
models, using the Kiel code JSteenbock & Holwegerl 11984 , 
with the model atom of iRentzsch-Holm (1996i) . The Kiel code 
uses a generalisation of the iDrawin I d 19691) formalism to take 
into account excitation and ionisation of the nitrogen atoms 
by inelastic collisions with neutral hydrogen atoms. A scaling 
factor, Sh, permits to reduce the efficiency of these collisions 
(Sh < 1), to switch them off (Sh - 0), or to consider them "in 
toto" (Sh = 1). Consis tent with the Kiel group and in line with 
iRentzsch-HoimI (Il996h , we favour Sh = 1/3, even though we 
computed NLTE corrections in addition for Sh = and Sh = 1 . 
In Fig. [3] the departure coefficients of the atomic levels relevant 
to the selected lines (see TablefTJ are shown for the (3D) model. 

The NLTE corrections obtained from the (3D) and the HM 
model are given in Table|4] They are added to the abundances 
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Fig. 3. Departure coefficients logZ?, - log(«NLTE/nLTE) of se- 
lected levels of neutral nitrogen as a function of optical depth 
logTsooo- For each line involved in our abundance analysis, the 
identification of the upper and lower level is listed in Table[T] 
The shown /?, refer to the (3D) model; similar results are ob- 
tained for the HM model. 



obtained from the LTE analysis based on the 3D and the HM 
model, respectively, to obtain the final NLTE nitrogen abun- 
dances. 



5.2. Influence of horizontal temperature fluctuations 

Using the (3D) model to compute departures from LTE is an 
extreme case in which all horizontal temperature fluctuations are 
neglected. To estimate how this approximation affects departures 
from LTE, we produced horizontal and time-averaged models by 
grouping columns into twelve bins according to their emergent 
continuum intensity in the vertical direction. A similar procedure 
was used by Aufd enberg et al. I (120051) to estimate the effects of 
horizontal, temperature inhomogeneities. 

These twelve components represent the different horizontal 
temperature structures associated with the presence of granules 
and intergranular lanes. Each component has a weight that de- 
pends on its surface area fraction. The cool downdrafts are as- 
sociated with the groups of lowest intensity, and the warm up- 
flows with those of the highest intensity. We computed NLTE 
corrections using each of these components, which were treated 
as standard plane-parallel model atmospheres (1.5D approxima- 
tion). The components allowed us to study the variations in the 
NLTE corrections as a function of the temperature structure of 
the flow where the line is formed. The results, for a selection of 
lines, are shown in Fig.|4] 

Lines with excitation energies higher than 11 eV exhibit lit- 
tle differences in departures from LTE for the twelve group- 
averaged models. For the lines with lower excitation energies 
(Flow ~ 10.3 eV), the groups corresponding to the cool down- 
drafts (groups 1 to 6) have a correction that is distinctly larger (in 
absolute value) than for the groups corresponding to the warm 
upflows (groups 7 to 12). We consider the variation in the NLTE 
corrections of the various components, relative to the result of 
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Fig. 4. For five representative lines the ID-NLTE corrections 
from the twelve group-averaged models, ordered according to 
increasing continuum intensity (radiation temperature) from left 
to right, are shown together with the result for the global (3D) 
model (rightmost, bold symbol). The lowest curve is plotted at 
the true ordinate level, while each of the others curves is shifted 
up by 0.05 dex with respect to the previous one for clarity. 
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Fig. 5. Results of the line profile fits for the 868.3 nm line. 
The synthetic spectra computed with SYNTHE and using the 
Holweger-Miiller model (thin lines) are superimposed on the ob- 
served spectra from the four solar atlases (KF - Kurucz flux, NF 
- Neckel Flux, DI = Delbouille Intensity, NI = Neckel Intensity; 
thick lines) 



the (3D) model, as an indication of the uncertainty associated 
with the ID approximation. 

The variations in the NLTE-corrections are small, even for 
the lower excitation lines, for which the difference between the 
extreme groups is less than 0.05 dex. We therefore expect that 
the results we obtained using the (3D) model will not disagree 
by more than 0.03 dex with the results of a full 3D-NLTE com- 
putation. 

6. Nitrogen abundance determinations 

6.1. LTE abundance from selected lines 

The following abundance analysis is based on a subsample of 
lines given in Table[Tl We determined individual LTE nitro- 
gen abundances, A(N)lte using the 3D, (3D), IDlhd, and HM 
model, respect i vely, and t he equivalent width m easurements by 
iGrevesse et alj (Il990h and lBiemont et al.l (1 19901) where available. 
The results are given in Table|2l In this Table and in the others 
following, we indicate the abundances to three decimal places, 
following the prescription of Bevington & Robinson ( 2003) of 
retaining one more significant digit than dictated by the error 
For nitrogen lines, which are formed in the deep photosphere, 
the 3D correction, defined as A(N)3d-A(N)idlhd' column (13), 
is small, in the range -0.024 - 0.008 dex. The "granulation 
abundance correction", A(N)3d-A(N)0d), (indicated in Col. (12) 
of Table|2]), which measures the influence of horizontal fluctua- 
tions, is somewhat larger, and negative for all the lines. We re- 
ported similar be haviour for the Oi 615.8 nm line in previous 
work (see lCaffau et al.ll2008h . 

6.2. Line-profile fitting 

We selected four clean and not too weak lines from the set of 
lines in Table[I] (746.8 nm, 821.6nm, 868.3 nm, and 1011.4nm) 
for which we fitted the line profile for all the four observed at- 
lases (two fluxes and two centre-disc intensities), using ID mod- 
els with SYNTHE as a line formation code. We discarded the 



line at 1 11 .4 nm for which we could not obtain a reliable fit; the 
abundance that we derived from this line was also too high to be 
consistent with those found for other lines, although we note that 
the four different spectra consistently indicated the same (high) 
abundance. The line list used in the profile fitting of this line 
appears to be incomplete. Nevertheless, we decide to retain this 
line in the abundance determination inferred from the EWs, be- 
cause there it is consistent with the other values, indicating that 
the line decomposition method used to determine EW is reliable. 

The A(N) derived from the line at 821.6 nm in the four dif- 
ferent spectra exhibits a scatter of 0.050 dex, when considering 
the (3D) -NLTE abundance, which is considerably larger than 
the scatter found for the other two lines. The difference between 
the A(N) derived from the two disc-integrated spectra is about 
0.06 dex, while the A(N) from the two centre-disc spectra differ 
by about 0.1 3 dex. Neither the two disc-integrated spectra nor 
the two centre-disc spectra do agree (see Fig.|2]l. The difference 
in EW is of the order of 15%. This is much larger than expected 
from the high S/N ratio of the spectra. 

The line fitting results are given in Table[3] An example of 
the excellent agreement achieved in the line fitting using the HM 
model, is shown in Fig.|5] 

The differences between the nitrogen abundance from both 
the fitting procedure and by matching the measured EWs can be 
explained by the influence of blending components, which are 
accounted for in different ways in the line-fitting approach and 
in the EW method, respectively. The different codes used may 
also be responsible for part of the difference. 

Since the majority of the nitrogen lines are blended, weak, or 
both, the line-fitting approach is, in principle, more appropriate. 
However, line fitting with 3D spectra can become a tedious task, 
and we therefore prefer to rely on equivalent width measure- 
ments. We believe that we cannot ach ieve superior measurem ent 
of the EWs than the c areful works by lGrevesse et alJ (1 19901) and 
iBiemont et al ] (119901) . In principle,' we could also measure the 
EWs in data from the other atlases, but our above line fitting 
exercise implies that this would introduce an additional scatter 



Table 2, A(N)lte from selected N i lines with EWs from the Uterature using LinforSD 
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Notes: Col.s with G are from lGrevesse et al] il990l) . with B are from lBiemont et al] ( Il990h . Col. (1) is the wavelength; col.s (2)-(3) the EWs; 
col.s (4)-(l 1) A(N) from 3D, <3D), IDlhd, and HM model, col.s (12)-(13) two different 3D corrections. 




Table 3. A(N) of N i 
SYNTHE 



from line fitting using ID models with 



746.80 746.82 



746.84 
(nm) 



746.86 746.88 



Fig. 6. The unblended 3D profile of Ni 746.8 nm is compared 
to the observed solar centre-disc profile of the Neckel Intensity 
atlas. 



of the order of at most approximately 0.03 dex, due to the differ- 
ences between atlases. As we shall see, this is much less than the 
scatter derived from the different lines. O ur final choice i s there - 
fore to derive A(N) from the EWs from Grevesse et al.l d 19901) 
and Biemont et al. ( 1990), using both 3D and ID models. 

As a consistency check, we also compared the 3D profiles 
with the observed spectra, as shown for one example in Fig.|6] 
The agreement is in general good, given that the blending line 
visible in the observed spectrum were not considered in the syn- 
thetic 3D profile. 

6.3. 3D-NLTE abundance 

Table|4] provides the NLTE corrections for both the (3D) and 
the HM models. Our final 3D-NLTE nitrogen abundance is ob- 
tained by averaging the 12 individual 3D-NLTE abundances 
(Col. (2)-i-Col. (4)), with equal weight. The individual as well as 
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Notes: Col. (1) is the observed spectrum identification, KF: Kurucz 
Flux, NF: Neckel Flux, NI: Neckel Intensity, and DI: Delbouille 
Intensity; col. (2) is the wavelength; col.s (3)-(5) the A(N)lte from line 
fitting with HM, ATLAS9, and <3D) model; col. (6) is the A(N)nlte 
for the (3D) model, obtained by adding the NLTE-correction computed 
with Sh = 1/3 to the abundances given in column (5). 



the averaged NLTE-corrections are given in Col.s (4), (5), and 
(6), respectively. The result is: 



A(N) = 7.85 + 0.12 
A(N) = 7.86 + 0.12 
A(N) = 7.87 + 0.12 



for Sh = 
for Sh =1/3 
for Sh = 1 



(1) 



With the EWs from iBiemont et~aLl (Il990l) . A(N) remains the 
same while the scatter is reduced to 0.06 dex. Since the scatter is 
insensitive to the choice of Sh due to these lines being formed at 
similar depths in the solar photosphere and having similar, and, 
in any case, small, NLTE corrections. 



Table 4. A(N)lte and AA(N)nlte from selected N i lines with disc-centre EWs from the literature obtained with LinforSD and Kiel 
code for NLTE computations. 
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average 7.885 7.890 -0.028 -0.020 -0.037 8.007 8.007 -0.032 -0.023 -0.041 



Notes: Cols with G are from lGrevesse et alj ( fT990h . with B are from lBiemont et al.l ( fT990h . Col. (1) is the wavelength; col.s (2)-(3) A(N) from 3D 
model, col.s (4)-(6) the corresponding NLTE corrections for Sh = 1/3, 1, 0, respectively; col.s (7)-(8) A(N) from HM model; col.s (9)-(ll) the 
corresponding NLTE corrections for Sh = 1/3, 1, 0, respectively. 



Table 5. Solar metallicity (Z and Z/X) for different choices of C 
and Ne abundance 



7.2. Estimating tine soiar metaiiicity 
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7. Discussion 

7.1. 3D effects and NLTE-corrections 

A decisive point in favour of using 3D hydrodynamical mod- 
els (or a semi-empirical model) is that there is no need to in- 
voke mixing-length theory. All the nitrogen lines used in our 
abundance determination arise from transitions between levels 
of high excitation energy, so that they are formed deep in the so- 
lar photosphere, in the range -1 < logr < +1, and hence are 
sensitive to the choice of the mixing-length parameter in classi- 
cal ID models. A change of ^ in the value of the mixing-length 
parameter of 1 .0 that we adopt, produces a change in the abun- 
dance with the IDlhd model of on average '^^^^^ dex. The "gran- 
ulation abundance correction", A(N)3d-A(N)(3d>, is always neg- 
ative, w -0.05 dex on average, in agreement with the results 
of iSteffen & Holwegej (1200 2). The total 3D correction, defined 
as A(N)3D-A(N)iDj^„Q(aMLT = 1-0) is only slighdy negative for 
most of the considered lines. Therefore, the 3D model provides 
a slightly smaller abundance (by « -0.01 dex) than the IDlhd 
model. 

The NLTE effects for the analysed lines are small and nega- 
tive, so that A(N)lte > A(N)nlte- The average NLTE correction, 
for Sh = 1 /3, is about -0.03 dex. 



Our favoured values for the solar nitrogen abundance, 
A(N)=7.86, and th e solar oxygen abundance, A(0)=8.76 
dCaffau et al.L l2008h . imply a revision of the solar metallicity 
Z. The two other elements important to the determination of the 
solar metallicity (in terms of fractional mass) are neon and car- 
bon. 

Neon is not measurable in the spectrum of the solar pho- 
tosphere. The usual procedure is to measure the ratio Ne/O 
in the solar corona, and assume this value to be the same as 
in ttie photosphere. We may consider Ne/O=0 .15 (As plund et 
al. ' 2005) or Ne/O=0.18 (Grevesse & Sauvall Il998h . imply- 
ing A(Ne)=7.94 and A(Ne)=8.02, respectively. By studying Ne 
abundances in nearby B-type stars iMorel & Butlen (12008.) found 
a mean A(Ne)=7.97 + 0.07, and suggest that this value should 
also be representative of the solar Ne abundance. 'Wang & Liul 
(2008) pointed out that the Ne/O ratios in Hn regions and plan- 
etary nebulae equal on average 0.25, and suggested that the cur- 
rent estimate of the solar Ne/O ratio could be too low. If the 
solar Ne/O ratio was so large, the Ne abundance would be 8.16. 
We may therefore consider four possible values of A(Ne), rang- 
ing from 7.94 to 8.16. For the carb on abundance, we ha v e thre e 
choices, A(C)^8 52 according to iGrevesse & Sauvall ( Il998l) . 
A(C)^8.39 from lAsplund et alTI (l2005|), or A(C)=8.592 from 
lHolwegeii(l2001h . 

The situation is summarised in Table|5] where our A(N) and 
A(0) are used with all possible combinations of A(Ne) and A(C) 
to provide Z and Z/X. The metallicity Z spans the range 0.0146 
to 0.0167, while Z/X spans the range 0.0199 to 0.0228. All these 
values are considerably higher than the values recommended by 
lAsplund et all (l2005l) (Z = 0.0122 and Z/X = 0.0165), and 
the increase is mainly driven by our higher oxygen abundance, 
although the adopted values of N, C, and Ne also play a role. 



7.3. Metallicity from helioseismology 

This upward revision alleviates the tension between photo- 
spheric abundances and helioseismic data, and the two may es- 
sentia lly agree. Using the solar relative co mposition of Asplund 
et al. ( l2005h and OP opacities. iBasu & A ntia (2008) inferred 
from the depth of the convection zone that Z/X = 0.0218 + 
0.0008, although the computation should be repeated with our 
solar composition, but we note that this value is already within 
the range of the Z/X in Table|5] Using the I Asplund et al. I (l2005h 
solar composition and one of their test models, IChaplin et al.l 
( 120071) derived Z = 0.0161 + 0.00008 from low-degree modes, 
again in the range of values of Table|5] Using the more exten- 
sive set of models o btained from a Monte Carlo simulation, 
IChapUn et all (|2007|) inferred values of between Z = 0.0187 
and Z = 0.0239. These values are definitely higher than any of 
the values in TablelSl From the solar sound-speed profile, Antia 
& Basu (2006) derive Z = 0.0172 + 0.002. Again, this value is 
higher than any value in Table|5] even though it is within 2cr of 
the highest value in the table. 



8. Conclusions 

Our recommended value of the solar nitrogen abundance is 
A(N)=7.86 + 0.12, which takes into account 3D effects and de- 
viations from LTE. This value is between the value inferred by 
[Asplund et al. (2005) (7.78 ±0.06) and that found by Grevesse & 
Sauval ( 1998) (7.92 + 0.06). Our A(N) is slightly lower, although 
consistent within errors, with the result found by Holweger 
(|20011A(N)=7.931 +0.111). 

In the course of our analysis, we found considerable differ- 
ences between the various solar spectral atlases, introducing ad- 
ditional uncertainties in the abundances derived from individual 
lines. The disagreement is by far larger than what can be at- 
tributed to the noise in the spectra. This fact has been pointed 
out already in our previous solar abundance studies. We believe 
that the quality of the available solar spectral atlases, which date 
back 30 or more years, should be improved by making use of 
up-to-date technology. 

Our most robust esti mate of the solar Z is obtain ed by adopt- 
ing the abundances of lGrevesse & Sauvall (Il998h for all ele- 
ments except for oxygen, nitrogen, and neon, for which we take 
our own abundances, A(0)=8.76, A(N)=7,86, and A(Ne)=8.02. 
The latte r results from the ratio Ne/O=0.18 given bv Grevesse 
& Sauval (Il998h . combined with our oxygen abundance. In this 
way, the resulting solar metallicity is Z = 0.0156 and Z/X - 
0.0213, respectively. 

Pending a new determination of the solar C abundance, we 
can say that our new results for N and O represent a significant 
step in forward reconciling photospheric abundances with helio- 
seismic constraints on Z. 
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